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Low-temperature magnetization process of the nanoscale iron cluster in linearly sweeped fields 
is investigated by a numerical analysis of time-dependent Schrodinger equation and the quantum 
master equation. We introduce an effective basis method extracting important states, by which 
we can obtain the magnetization process effectively. We investigate the structure of the field 
derivative of the magnetization. We find out that the antisymmetric interaction determined from 
the lattice structure reproduces well the experimental results of the iron magnets and that this 
interaction plays an important role in the iron cluster. Deviations from the adiabatic process are 
also studied. In the fast sweeping case, our calculations indicate that the nonadiabatic transition 
dominantly occurs at the level crossing for the lowest field. In slow sweeping case, due to the 
influence of the thermal environment to the spin system, the field derivative of the magnetization 
shows an asymmetric behavior, the magnetic Fohn effect, which explains the substructure of the 
experimental results in the pulsed field. 

KEYWORDS: Nanoscale magnet, magnetization process, time-dependent Schrodinger equation, quantum master 
equation, Lanczos diagonalization, Runge-Kutta method, nonadiabatic transition 



i §1. Introduction 



Nanoscale magnets have attracted much interest from 
the point of view of their quantum dynamics.^ As such 
magnetic materials, Pe„ (n — 6, 10, 12, 18) have been 
synthesized. In these compounds, each iron ion Fe 3+ 
plays a role of an S = 5/2 spin and forms a ring- 
structure molecule. Each of these compounds reveals 
a step-like magnetization process M{H) at low temper- 
atures. These steps appear due to the fact that the 
energy levels are discrete, which is one aspect of the 
quantum nature. First, let our attention be focused on 
the magnetizations differentiated by the magnetic field, 
dM/dHa> At magnetic fields where the magnetization 
grows rapidly, dM/dH shows peaks with various heights. 
In two of J-hc compounds [M Fe 6 (OCH 3 )ia.(dbm) 6 ] + 
(M=Na,Li)0) and [Fe(OMe) 2 (O 2 CCH 2 Cl)]i J§ , which 
are abbreviated by Fe6 and Feio respectively, the peak 
at the lowest magnetic field has the largest height; the 
heights of the peaks get smaller gradually at higher fields. 
How does the difference of the peak heights appear? 

The main structure of the magnetization process is 
determined by the Heisenberg-like exchange interactions 
between the iron ions, where the z-component of the 
magnetization M z is a good quantum number. Eigcn- 
states of various values of M z show many crossings in 
the energy-level structure as function of the magnetic 
field H. The ground state is given by the lowest eigen- 
state at the given field, which changes piecewise. The 
typical case of antiferromangetc interactions is shown 
schematically in Fig. ??. The singlet ground state, which 
is nonmagnetic, is realized at H = 0. At H = Hi and 
H2, the ground state changes, where the magnetization 
of the ground state increases. If a system includes only 
the Heisenberg-like exchange interactions, dM/dH re- 



veals the d function at the crossing. In experiments, 
however, they have finite peaks. There are two possi- 
ble origins for the finite heights. One is an effect due 
to finite temperature. The other is an effect due to 
some quantum-mechanical operators. Let us consider 
the case of the thermal effect. At low temperatures, 
only the two eigenstates concerning the crossing are rel- 
evant. Main contribution in the partition function con- 
sists of the Boltzmann weights of these states, and the 
thermal behavior is characterized by the energy differ- 
ence between them. This effect causes almost equiva- 
lent peaks at all the crossings in dM/dH, which is not 
the case in the experiments. Therefore, it seems that 
the origin of the difference between the peak heights is 
due to some operators which do not commute with M z . 
What kind of operator does induce such a transition? 
In this paper, we investigate such interactions. We con- 
sider the lattice structure of the iron clusters and find 
the existence of the antisymmetric interaction, namely 
the Dzyaloshinsky-Moriya interaction. This interaction 
which does not commute with M z is found to reproduce 
well the feature observed experimentally in the magneti- 
zation process of the iron clusters. 

It is known that when the parameter is swept cross- 
ing the transition .RQip.t, the system shows the nona- 
diabatic transitions ErEP The transition rate depends 
on the sweeping speed of the field. In particular, the 
present case belongs to the so-called Landau-Zener- 
Stiikclberg (LZS) mechanism. Tunneling phenomena of 
the metastable magnetization of the Ising-type model 
were-studied in the viewpoint of the nonadiabatic transi- 
tion© Here we treat a model for the nanoscale magnets 
and clarify how effects of the nonadiabatic transition ap- 
pear. 

The peaks in dM/dH of Feg and Feio are known to 
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be asymmetric in higher- and lower-field regions with re- 
spect to the center of each peak. In the receding side of 
the pulsed field, the slope of the peak in dM/dH is gen- 
tle in comparison with that of the approaching side. In 
the extreme case, a satellite peak is accompanied in the 
receding side. Similar behavior of hysteresis in the mag- 
netization process is known in [Fe(salen)Cl]2 (s£tcalled 
Fe 2 )B> and K 6 [Vi5As 6 042(H 2 0)] (so-called V W )S3 The 
behavior of the vanadium compound was investigated 
as influence of the thermal environment which is called 
magnetic Fohn effectO In this work, we investigate the 
magnetization process of Fe6 with the effect of environ- 
ment using the quantum master equation. We find out 
that the asymmetry of the peaks of the iron clusters can 
be understood as an effect from the environment. 

This paper is organized as follows. The next section 
is devoted to the model Hamiltonian and the calculation 
method. The antisymmetric interaction in the model 
is introduced, which explains the main structure of the 
peaks in dM/dH. As to the method, we introduce an 
effective basis method in a dynamical approach in the 
Schrodinger equation and also in the quantum master 
equation. In §3 the magnetization process at zero tem- 
perature is discussed. In §4 the magnetization process at 
finite temperatures is investigated. Concluding remarks 
are given in §5. 

§2. Hamiltonian and Method 

2.1 Spin Hamiltonian 

We consider a system that has a large interaction of 
Heisenberg type between two neighboring spins. In addi- 
tion, there presumably exists other anisotropic interac- 
tions determined from the lattice structure of the com- 
pound even though they are smaller than the Heisenberg- 
type one. Let us consider the lattice structure of the iron 
clusters. In both Fe6 and Feio, two neighboring iron ions 
have a mirror symmetry with respect to the plane per- 
pendicular to the line connecting these iron ions. In this 
case, it is known that the Dzyaloshinsky-Moriya (DM) 
interactions 

N B 

HuM = Y^D j -[S j xS i+1 ], (2.1) 

3=1 

may be present, and that Dj is parallel to the mirror 
plane. Here, N s denotes the number of sites and Sj de- 
notes an S = 5/2 spin operator at site j. The system is 
periodic, i.e., Sn s +i = Si. It is known that, in every two 
neighboring pairs, branches of atomic bonds are alternat- 
ing. Therefore, Dj and Dj+i have opposite directions. 
In this work, we treat the regular hexagon (N s = 6) with 
alternating Dj 

D j = (-iy(D,0,0), (2.2) 

as shown in Fig. ??. Thus, as a candidate model for the 
system of Fe6 and Feio, we examine the Hamiltonian 

H sys = Tto + Hbm, (2.3) 
H =Y t [JS j -S j+1 -h(t)S<], (2.4) 

3=1 



where J is an amplitude of antiferromagnetic exchange 
interaction, and h(t) is the external magnetic field. Here 
we consider the case of linear sweeping, h(t) — ct. In 
this work we examine the case where D is much smaller 
than | J\. A main structure of energy is determined by 
Ho as a function of the field, as shown schematically in 
Fig. ??. Because J ^> D, the effect of Hbm appears only 
in the vicinity of a level crossing, and causes the avoided 
level crossing structure. 

2.2 Calculation method 

Let us obtain the magnetization process. In this pa- 
per, we introduce a dynamical approach to obtain the 
magnetization process. It is possible in principle to cal- 
culate the adiabatic magnetization process M(h) of H sys 
by diagonalization. However, in order to obtain smooth 
curves, very fine steps in the field are necessary, which re- 
quires a large number of diagonalizations. Moreover this 
dynamical approach provides information on the dynam- 
ical processes which cannot be obtained by diagonaliza- 
tion. 

To obtain the magnetization process in the ground 
state (T = 0), we solve the time-dependent Schrodinger 
equation 

i7L^|V)=W sys |V>>, (2-5) 

numerically. This equation is made dimensionless by the 
scaling t — Jt/h. We take J as an energy unit; h = 
h/J and D = D/J are independent parameters in the 
Schrodinger equation. The external magnetic field h is 
increased linearly from zero, which corresponds to the 
pulsed field in experiments. In particular, although the 
pulse velocity has a sinusoidal shape in experiments, here 
we approximate it by a constant velocity. As far as the 
velocity is small, we can obtain the adiabatic process in 
both cases. The antiferromagnetic ground state for h = 
is taken as an initial state. 

To treat the thermal environment, we have to consider 
the effect of environment, where the total Hamiltonian 
is given by 

Htot — H S ys + Hint + Hbath, (2.6) 

where 7ibath denotes a Hamiltonian of the heat bath. 
Hint denotes a term of interaction, which plays a role of 
interface exchanging between the heat bath and the spin 
system. We take the heat bath as the free boson system, 
namely, H ha .th = Z) a ^ w a( & l & " + 1/2). H m t is taken 
as Hint = ^V / ^S Q (^J l + b a )X, where X is an operator 
in the spin system. Here, {4>m\X\^)m') is taken as a 
unit when \M — M'\ = 1 and vanishes otherwise. Under 
the condition of small coupling A, the quantum master 
equation for the density matrix 

^ = ^[7W(t)]-rp(t), (2.7) 

has been derived.© Here, 

Tp{t) = {[X, Rp(t)] + [X, Rp{t)]t}\ 2 K/h, (2.8) 
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(k\R\n) = \xk,n<5> (^-jr 



$(w) 



(2.9) 



(2.10) 



where (3 is the inverse of temperature, namely k^T — 
1/(3. Here, \k) and \n) represent eigenstates for H sys 
with energy eigenvalues and E n , respectively. If we 
assume that a main contribution of the heat bath comes 
from the lattice, the bosons in the environment are con- 
sidered to be phonons. Then, the spectral density I(uj) 
is given by Iqlo 2 9{uj) where 9(uj) is the step function. We 
solve the quantum master equation (2.7) to obtain the 
magnetization process at finite temperatures. 



2.3 Effective basis method 

Since the number of states A D = (25* + 1)^=46656 
is large, solving the Schrodinger equation (^) using all 
the states in the system is difficult, in particular for long 
t ime s. The calculation of the quantum master equation 
(2.7) is even more difficult because the number of el- 
ements of the density matrix is N^. Therefore, some 
approximation method which extracts important states 
of the system is required. 

Now we introduce an approximate method for this pur- 
pose. First note here that many eigenstates are irrelevant 
to study the magnetization growth with the increasing 
field in the ground state through the quantum transitions 
of the present interests. In spite of this fact, the irrele- 
vant states and important states are treated equivalcntly 
in full-size calculations, and this is inefficient. Now let us 
consider a method to extract the important states which 
contribute to the quantum transitions. 

Here, the following procedure is performed, which 
we call an effective basis method. At first, one ob- 



tains the ground state |-0g) of TC S y S for a field h(t) at 
a time t (initially, t — 0) by the Lanczos diagonalization 
in the basis taking the z direction as the quantization 
axis. Note that for the present system with dimension 
Njj = 46656, this is still possible. From IV'g), one ex- 
tracts \<f>M.) = \<t>M M )/y {<t>M,\4>M g ), where 



(2.11) 



where |£fc) denotes the fc-th basis state with the magneti- 
zation M z , i.e., Sl ot \£k) = M*|&). Here, 5*°* is defined 
by J2f=i Sj' Note that M z is an integer ranging be- 
tween zero and the saturation value, which is equal to 
15 (=| x 6) for the present case. The new state \<j>M z ) 
is a representative of the states of M z . We have checked 
that the set of states {|0m 2 )} represents well the behav- 
ior of Hsys at low energy. Using this set of states, one 
constructs a small matrix 



ttsys — 



(0o|W s y S |0u) 

(</>l|W S y S |0u) 

(015|?4ys|<A)) 



&)|W s ys|0l) 

1 Ways 1 01/ 
^l^sysl^l 



^o|W sys |0i 5 ) 

h | H SyS | 015 ) 



15 |^sys|015) 



(2.12) 

Here, we take only the ground state to generate the effec- 
tive basis. In this paper, we find that even in the present 
basis, we obtain good results. However, we can improve 
the approximation by adding more basis states extracted 
from a few low-energy states. 

In this work, we perform Rungc-Kutta calculations for 
the small matrix 7i sy s instead of the full-size matrix 7i S ys- 
In the Runge-Kutta calculations, the normalization of 
the wave function deviates easily from unity. We have 
checked it carefully and used a small enough time step. 
The initial condition is |0o)- To keep the quality of the 
approximation the same, we upd ate a set {|0a/. )} by the 
above mentioned method, i.e. ( 2.11 ), after time evolu- 
tion by Ti. for a while. It is found that {|0m z )} gives 
a good approximation if it is updated as often as the 
level crosses. We have verified that no significant differ- 
ences appear when the frequency of the update is varied. 
When the sweeping speed is small and the system be- 
haves adiabatically, the states with the lowest and the 
second lowest energies are important. The method de- 
scribed here gives the magnetization process M(h) as the 
adiabatic process very effectively. Actually, this dynam- 
ical method is a natural way corresponding to experi- 
ments of pulsed magnetic field. 

As we mentioned, the present method provides infor- 
mation of nonadiabatic process for fast sweepings, that 
cannot be obtained by a static Lanczos diagonalization 
only. To demonstrate the validity of the above approx- 
imation for the nonadiabatic transition, we have per- 
formed the calculation of the Schrodinger equation for 
S = 1 for D = 0.01. In this case, the first avoided level 
crossing has an energy gap AE/J ~ 0.0196. According 
to the LZS mechanism, the transition rate at an avoided 
crossing is given by 

n(AE/J) 2 



P = 1 



exp 



2\M Z - M' z \ 

"-max/ f-max 



(2.13) 



when the magnetic field is linearly sweeped from zero 
to /imax during a time of i max . For the above energy 
gap, h max = 6, ? max = 3000, M z = 1 and M' z = 0, the 



transition rate P given by (2.13) is 0.26. After the corre- 



sponding level crossing, the present approximation gives 
for this rate 0.26. This agreement indicates the valid- 
ity of the above approximation. For the quantum mas- 
ter equation, we use the same effective basis. We have 
also checked the agreement of a result from the quantum 
master equation for a full-size Hamiltonian with a small 
dimension and that from the quantum master equation 
for the Hamiltonian with its effective basis, which also 
indicates the validity of the approximation in the quan- 
tum master equation. 

§3. Magnetization Process at T = 

Figure ??(a) shows the result of the magnetization 
process for f max = 30000, h max — 15 and D — 0.01 
together with the complete stair case of the ground state 
of the model of D = 0. Every step indicates the position 
of a level crossing at which the magnetization grows from 
M to M + 1. Because a final value of the magnetization 
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reaches full saturation, the process is almost adiabatic at 
every crossing. Figure ??(b) shows dM/dh for D = 0.01. 
Peaks with various heights induced by the DM interac- 
tion are formed at the level-crossing positions denoted 
by crosses. It is noticable that the heights of the peaks 
are different from each other. Especially, the first peak 
is the highest. The heights gradually become small as 
the field is increased although the heights in a few last 
peaks get higher again. These features of the different 
peak heights agpee well with, the experimental results in 
dM/dH for Fe<0 and Feio iff 

Single-ion-type anisotropy and dipole-dipole interac- 
tions may causes broadening of dM/dH . However, the 
contribution of the single-ion anisotropy is considered to 
be small because the orbital angular momentum vanishes 
in the each iron ion Fe 3+ of the above materials. The off- 



diagonal element 
interactions, 



&2M+i|%lip|02Af) of the dipole-dipole 



i,3 



(3Mb) 1 



r Si • Si 



OOSi-ry)], (3.1) 



vanishes because \4>2hi) and \4>2M+i) have wave numbers 
7r and 0, respectively. This means that contribution of 
the dipole-dipole interactions is also small. Thus, we 
conclude that the transition due to the alternating DM 
interactions is the major origin of the peaks of dM/dH 
of Fc6 and Feio with different heights. 

On the contrary the dM/dH curve of-Fei2, which is 
an abbreviation of [Fe(OCH3) 2 (dbm)]i 2 f$ reveals that 
the height of the first peak in dM/dH is not the largest 
and that the second one is larger than the first one.li-^ It 
was reported that the lattice structure of this material 
has a large deviation from the regular polygon, while in 
Fe6 and Feio, the lattices are almost regular polygons. 
In order to take the deviation into account, we study 
a lattice with locally non-uniform magnetic interaction 
{Ji}- In our calculations, however, the first peak is still 
the highest in the system with various sets of { Ji}. 

In the dimer material Fe2, the second peak is higher 
than the first one.B The interactions in these compounds 
should be clarified in the future. 

So far we studied the adiabatic behavior for slow 
sweeps, which probes the magnetization processes in the 
ground state. Next, let us consider the dynamical ef- 
fect of the sweeping of the magnetic field in the case of 
a larger speed. The result for i max = 7500 is shown in 
Fig. ??. The magnetization at a field where the satura- 
tion value is reached in Fig. ??(a) becomes smaller. This 
means that there are some level crossings at which nona- 
diabatic transitions occur. At h = 0.7, the magnetization 
jumps up to M = 1 in the adiabatic process. However, 
in Fig. ??(a), the step height is clearly lower, suggesting 
a nonadiabatic transition at the first level crossing. In 
the fields after the first crossing before the second one, 
the state is given approximately by a linear combina- 
tion of \<fro) and \<f>x), i.e., = c |</>o) +ci\(f>i), where 
\<f>o) and \(f>i) behave almost as eigenstates for the field 
away from the crossings. The probabilities of \<po) and 
\<t>i) are found to be |c | 2 = 0.22 and | d | 2 = 0.78, re- 
spectively. The plateau height m satisfies the relation 



m = \ci\ 2 = 1 — | co | 2 = 0.78. For fields between the sec- 
ond and the higher-field crossings, the growth of the mag- 
netization is briefly linear. The population of | cx | 2 = 0.78 
stays in the ground state, which gives the linear increase 
of the magnetization. The saturated magnetization for 
t max = 7500 agrees with 15 xm. These facts indicate 
that the system shows almost adiabatic process at cross- 
ings other than the first one within this sweeping speed 
and that the energy gap at the first crossing is found to 
be the smallest. Thus, the transition at the first cross- 
ing plays a role of a bottleneck process for the growth 
of the magnetization. Note here that the population at 
\4>o) does not cause changes in the magnetization process 
within the present approximation although it is also ex- 
posed to nonadiabatic transitions in higher fields in the 
full-size Hamiltonian. To obtain such nonadiabatic tran- 
sitions, an improvement of the method of calculation is 
required. 

If the nonadiabatic transitions are experimentally ob- 
served as a quantum effect, we could capture the phe- 
nomena of the states which branch at the level crossings. 
Duration times i max = 7500 and i max = 30000 corre- 
spond to 3~6 [ns] and 11^22 [ns] in real time, respec- 
tively, by using J/fce = 10 ~ 20 [K] estimated from the 
susceptibility fit of Fe6 and Feio- In the present experi- 
ments, on the other hand, the duration time of the pulsed 
field is in order of 10 [ms]. In order to find the nonadia- 
batic phenomena in Fee and Feio, either a faster sweep- 
ing speed or a value of D smaller than 0.01 is necessary. 
When D is small, it is expected that D is proportional 
to the energy gap at the level crossing. The amplitude of 
the gap determines whether the transition is adiabatic or 
nonadiabatic. Then, for the appearance of the nonadia- 
batic features D/k-Q must be in order of 0.1 [mK], if the 
duration time is 10 [ms] and J/k-Q = 10 [K]. An exper- 
imental estimate of the amplitude of the antisymmetric 
interactions in Fe6 and Feio would be very helpful. 

§4. Magnetization Process at Finite Tempera- 
tures 

In the previous subsection, we studied the dynamical 
magnetization process in the pure quantum-mechanical 
way. At finite temperatures, however, influence from the 
surrounding environment may have significant effects. 
Here, we study the effect of the thermal environment 
on the iron clusters. 

Figure ?? shows the result of the magnetization pro- 
cess for the case of D = 0.01, A = 0.005, h max = 12 
and i max = 60000 in the condition of temperature 
ksT/J = 0.05. A satellite peak in dM/dh, is clearly 
observed around h = 0.9 in the receding side of the field 
near the main peak at h = 0.7. This satellite peak cor- 
responds to a dent of the magnetization process at an 
edge of the magnetization step. This satellite peak is 
considered to be a consequence of the magnetic Fohn ef- 
fect The key mechanism of the magnetic Fohn effect 
is the difference between two kinds of speed characteriz- 
ing properties of the total system. One is a speed deter- 
mined by the LZS mechanism, below which the quantum 
transition is almost adiabatic. The other is a speed con- 
trolled by Hint, above which the supply of energy from 



Magnetization Process of Nanoscale Iron Cluster 



5 



the heat bath is not fast enough to keep the system in the 
equilibrium during the change of the field. If the sweep- 
ing speed of the field is between these two, the magnetic 
Fohn effect occurs. This effect has been proposed as a 
general phenomenon in the processes which are almost 
adiabatic but with a little inflow of heat. til Because the 
process is almost adiabatic, the system is almost in the 
ground state, that is in the lower-energy state. How- 
ever, due to the interaction with the environments, some 
population is excited to the higher-energy level. But it 
is not enough to realize the thermal equilibrium distri- 
bution. This excitation occurs as far as the equilibrium 
population at the higher-energy level p%q(<x e~ ,3A (' 1 )) is 
larger than the dynamical one of the system p( 2 '(t), when 
A(/i) is the energy difference between the states. When 
the field gets away from the crossing point, the energy 
difference increases; pi^ decreases rapidly and becomes 
smaller than p^ 2 \t). There, the population at higher- 
energy level begins to relax to the lower-energy level. In 
this way, p^ (t) has a maximum at a field h in the reced- 
ing side from the crossing point. The behavior of p^ (t) 
causes a dent of the magnetization process. 

In the receding side of the field near the main peak at 
h = 1.4, a smaller swell appears around h = 1.6. In this 
case, the satellite due to the thermal redistribution is 
being merged into the main peak. Situations of merging 
depend on the distance between the satellite and the cen- 
ter of the main peak, and on the width of the main peak. 
The former is determined by the details of the heat bath 
and those of the interface interaction. The latter is given 
quantum-mechanically by the energy gap at the avoided 
level crossing. In the present case, the merging behavior 
is induced because the width of this main peak becomes 
larger in higher fields. This effect causes a gentle slope in 
the receding side compared with that in the approaching 
side. Thus, the asymmetry of the experimental peak in 
the field derivative of the magnetization is understood 
as a consequence of the magnetic Fohn effect induced by 
the thermal environment. 

Let us consider the temperature dependence of the 
satellite peak. Results are shown in Fig. ??. As the 
temperature is decreased, the satellite peak gets smaller 
and almost vanishes at k^T j J = 0.01. It is character- 
istic that the satellite peak disappears without merging 
with the main peak. This temperature corresponds to 
0.09 ~ 0.18 K if we take a experimental estimation of 
J/k-Q to be 10 ~ 20 K. This means that in the exper- 
iments for Fe6 at such low temperatures, the satellite 
peaks in dM/dH would disappear in the way we ob- 
served. A detailed comparison with such an experiment 
is required. 

§5. Concluding Remarks 

We have studied properties and the origin of the quan- 
tum transitions observed in the magnetization process of 
the nanoscale iron clusters at T = and at finite tem- 
peratures. We have proposed an approximate procedure 
to make the numerical analysis possible and obtain the 
magnetization process for the system where the dimen- 
sion of the basis states is large and conventional method 



is hardly applicable. We have found that the antisym- 
metric interaction reproduces the characteristics of the 
heights of the peaks in dM/dH at T = 0. The effect of 
nonadiabatic magnetization process was also discussed. 
The first level crossing plays a role of the major branch- 
ing point. At finite temperatures, the asymmetric be- 
havior of the peaks in dM/dH is induced by the thermal 
environment influencing the spin system as the magnetic 
Fohn effect. In experiments, observation of the mag- 
netization process is a powerful and direct method to 
estimate the coupling amplitude J of the Heisenberg- 
like exchange interactions. However, the magnetization 
process yeilds not only information on J but also other 
information characterizing the properties of the materi- 
als. In this paper, we have clarified such aspects of the 
iron clusters. 
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